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ABSTRACT 

It has been proposed that a galaxy’s nova rate might be enhanced by the production of nova progenitor 
binaries in the dense cores of its globular clusters (GCs). To explore this idea, relative nova rates 
in three Virgo elliptical galaxies, M87, M49 and M84, which have significantly different GC specific 
frequencies {Sn) of 14, 3.6, and 1.6, respectively, were measured over the course of 4 epochs spanning 
a period of 14 months. To simplify the analysis, observations of the nearly equidistant galaxies were 
made on the same nights, with the same integration times, and through the same filter (Ha), so that 
the relative numbers of novae discovered would reflect the relative nova rates. At the conclusion of 
our survey we found a total of 27 novae associated with M87, 37 with M49, and 19 with M84. After 
correcting for survey completeness, we found annual nova rates of 154^^g, I 89 I 221 ^od 95^]^4, for 
M87, M49, and M84, respectively, corresponding to ib-band luminosity-specific nova rates of 3.8± 1.0, 

3.4 ± 0.6, and 3.0 ± 0.6 novae per year per 10^° Lk,q- The overall results of our study suggest that 
a galaxy’s nova rate simply scales with its luminosity, and is insensitive to its GC specific frequency. 

Two novae, one in M87 and one in M84, were found to be spatially coincident with known GCs. After 
correcting for the mass fraction in GCs, we estimate that novae are likely enhanced relative to the 
field by at least an order of magnitude in the GC systems of luminous Virgo ellipticals. 

Keywords: novae, cataclysmic variables, globular clusters — galaxies: individual (M87, M49, M84) 


1. INTRODUCTION 

Classical novae are a subclass of cataclysmic variables 
in which a progenitor system increases in brightness by 
~I0-20 mag over a period of days to weeks and then re¬ 
turns to its pre-outburst magnitude over a period of days 
to years. The progenitor system is a short period, semi¬ 
detached binary consisting of a white dwarf primary ac¬ 
creting material from a Roche lobe filling late-type com¬ 
panion. If the accretion rate is sufficiently low, the ac¬ 
creted material forms a thin shell of degenerate matter 
on the surface of the white dwarf. Once the temperature 
and pressure at the base of the accreted envelope be¬ 
come sufficiently high, a thermonuclear runaway ensues, 
resulting in a nova eruption. Novae can reach an abso¬ 
lute magnitude, My ~ ~10, making them among the 
most luminous explosions in the Universe, easily visible 
in galaxies as distant as the Virgo cluster (see Shatter 
2008, for a review). 

The nature of the progenitor system suggests that no¬ 
vae will be recurring events. Once a nova eruption occurs 
the material will begin to build up again and the process 
will repeat. Recurrence times can be as short as I year 
(e.g., see Darnley et al. 2014; Henze et al. 2014) all the 
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way up to 10® yr, or perhaps longer. Despite the fact 
that all novae are believed to be recurrent, only novae 
with more than one recorded outburst are referred to as 
“recurrent novae”. Over time, the orbital period of a 
nova progenitor binary decreases as a result of angular 
momentum losses due to magnetic stellar winds from the 
tidally-locked donor star, or from gravitational wave ra¬ 
diation, or both. As the orbital period decreases, so does 
the average accretion rate onto the white dwarf primary, 
resulting in an increase in the recurrence time (Iben et 
al. 1992 ). In addition, in a given stellar population, the 
average mass of the primaries in nova systems are ex¬ 
pected to decrease with increasing time elapsed since the 
zero-age main sequence population formed (Tutukov & 
Yungelson 1995). A less massive primary needs to accrete 
more material to reach the necessary temperature near 
its surface to undergo a thermonuclear runaway (e.g., see 
Townsley & Bildsten 2005). Such systems will undergo 
less frequent (and dimmer) nova outbursts compared to 
systems with higher mass primaries (Ritter et al. 1991; 
Livio 1992; Kolb 1995). 

Accounting for the above factors, population synthesis 
models have predicted that K-band luminosity-specific 
nova rates, vki should be significantly higher in the 
younger populations found in spiral disks compared with 
older populations found in elliptical galaxies and in the 
bulges of spiral galaxies (e.g., Yungelson et al. 1997). De¬ 
spite these predictions, many nova surveys have found 
that, within observational uncertainties (which can be 
large), vk appears to be relatively independent of Hub¬ 
ble type (e.g. Ciardullo et al. 1990a; Ferrarese et al. 2003; 
Williams & Shatter 2004; Coelho et al. 2008; Giith et al. 
2010; Franck et al. 2012), although some surveys have 
found evidence that uk varies between galaxies, though 
not necessarily as predicted by models (e.g. della Valle 
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et al. 1994; Shara & Zurek 2002; Madrid et al. 2007). 
In addition, nova surveys of M31, in which nearly 1000 
novae have been observed (e.g., see Pietsch et al. 2007), 
have repeatedly found that the nova spatial distribution 
is more centrally concentrated than the background light, 
indicating a higher vj^ in the bulge than in the disk (e.g. 
Ciardullo et al. 1987; Capaccioli et al. 1989; Shatter & 
Irby 2001; Darnley et al. 2004, 2006). More recently Neill 
& Shara (2004) have found a similar result for M81. In 
light of the apparent discrepancy between theory and 
observation found in their early M31 survey Ciardullo et 
al. (1987) suggested that a significant number of nova 
progenitor binaries may have formed in M31’s globular 
cluster (GC) population and been subsequently injected 
into the bulge through 3-body encounters within clusters, 
or through tidal disruptions of entire clusters, or both, 
thereby enhancing the bulge nova rate. 

It has long been known that the number of Galactic 
low-mass X-ray binaries per unit mass are enhanced in 
GC environments (Clark 1975; Katz 1975), and it has 
been suggested that many Galactic bulge X-ray sources 
could have had their origin in GGs (Grindlay 1988). 
A similar enhancement of X-ray sources has also been 
seen in M31’s GG population (Grampton et al. 1984; 
Di Stefano et al. 2002). In recent years, it has been 
suggested that a significant population of close binaries 
with white dwarf accretors (e.g., nova progenitors) could 
also be formed in the dense cores of GGs (e.g., Pooley 
& Hut 2006). Observations by the Hubble Space Tele¬ 
scope (HST) and Chandra had already begun to reveal 
these binaries in Galactic GGs (e.g., Edmonds et al. 2003; 
Heinke et al. 2003; Knigge et al. 2002; Pooley et al. 2002), 
although there is no compelling evidence that they are 
enhanced relative to the field. 

Additional evidence for the idea that nova progenitors 
could be produced in GGs has come from the study of 
nova rates in Virgo elliptical galaxies. Shara & Zurek 
(2002) reported that the nova rate in the giant elliptical 
galaxy M87 could be as high as ~ 300 yr“^, while for 
the slightly more luminous elliptical, M49, Ferrarese et 
al. (2003) found a rate of i? ss 100. It is interesting that, 
like the ratio of observed nova rates, M87 has a specific 
frequency of GGs about three times that of M49 (Brodie 
& Strader 2006). Although, these results are consistent 
with the possibility that the large GG population of M87 
enhances its nova rate, absolute nova rates extrapolated 
from observed nova rates are largely uncertain given the 
small spatial coverage of these HST surveys, and compar¬ 
isons of these rates between surveys can be unreliable. 

In 2011 we began a definitive test of the putative rela¬ 
tionship between a galaxy’s GG specific frequency, S'^r, 
and its luminosity-specific nova rate, i^k- Specifically, we 
searched for novae in three essentially equidistant Virgo 
elliptical galaxies having different GG specific frequen¬ 
cies. The galaxies were observed under nearly identical 
conditions using the same telescope, observing cadence, 
filters, and exposure times. With this approach we were 
able to simply compare the relative numbers of novae dis¬ 
covered in each of the three galaxies directly to estimate 
relative nova rates, while avoiding the large uncertainties 
that result when attempting to measure absolute nova 
rates for different galaxies in different surveys. Prelim¬ 
inary results of our program were discussed in Shatter 


Table 1 

Galaxy Properties 


Galaxy 

Sn 

(m - M) 
(mag) 

(mag) 

(B - V)o 

(mag) 

{V-K)o 

(mag) 

M87 

14 

31.03 

9.49 

0.96 

3.21 

M49 

3.6 

31.06 

9.28 

0.96 

3.34 

M84 

1.6 

31.32 

9.99 

0.87 

3.28 


Note. — Sn are from Brodie & Strader (2006), (m — M) 
are from Tonry et al. (2001), Bt,o and (B — f/)o are from de 
Vaucouleurs et al. (1991), and {V — K)o are from Frogel et 
al. (1978). 


Table 2 
Observations 


Epoch 

Date 

Galaxy 

Exp 

Seeing 

Weather 

1 

2011 Feb 02 

M49 

12 X 5m 

0.99" 

phot 


2011 Feb 02 

M84 

12 X 5m 

0.82" 

phot 


2011 Jan 30 

M87 

12 X 5m 

0.87" 

phot 

2 

2011 Mar 30 

M49 

12 X 5m 

0.82" 

phot 


2011 Apr 02 

M84 

12 X 5m 

0.89" 

phot 


2011 Mar 29 

M87 

16 X 5m 

0.95" 

P cldy 

3 

2012 Feb 26 

M49 

11 X 5m 

1.19" 

phot 


2012 Feb 27 

M84 

12 X 5m 

0.72" 

P cldy 


2012 Feb 26 

M87 

12 X 5m 

0.96" 

phot 

4 

2012 Mar 21 

M49 

15 X 5m 

0.72" 

P cldy 


2012 Mar 22 

M84 

12 X 5m 

0.66" 

P cldy 


2012 Mar 21 

M87 

13 X 5m 

0.81" 

P cldy 


et al. (2013) and Gurtin et al. (2014). Here, we present 
the final results of the first phase of our ongoing Virgo 
galaxy nova study. 

2 . OBSERVATIONS 

We chose three Virgo ellipticals for our study: M87 
(NGG 4486), M49 (NGG 4472) and M84 (NGG 4374). 
As members of the Virgo cluster, these three galaxies 
are approximately equidistant. They also have similar 
luminosities, but they have widely varying GG specific 
frequencies, Sn (see Table 1). We observed all three 
galaxies over four epochs from February 2011 to March 
2012 using the Canada-France-Hawaii 3.6-m telescope 
(GFHT). All observations were made using the Mega- 
Gam which consists of a mosaic of 36 2048x4612 pixel 
GGDs covering a full field of view of ~ 1 deg^ at a res¬ 
olution of 0.187"/pixel (Boulade et al. 2003). Exposure 
times were kept nearly constant between galaxies, vary¬ 
ing only slightly due to weather in a small number of 
cases to ensure that the derived relative nova rates would 
be directly comparable. The observations are summa¬ 
rized in Table 2. 

We chose to image all three galaxies through a 
narrow-band Ha filter redshifted to the mean velocity 
of the Virgo cluster (Ac = 6584A, FWHM = 76A[~ 
3500 km s“^]). For nova surveys imaging in Ha has two 
principal advantages over broad-band imaging. Novae 
exhibit strong Ha emission lines shortly after outburst 
that fade slowly relative to the continuum, typically re¬ 
quiring months to decline by more than 2 mag. The slow 
decay rate in Ha means that we can detect novae over 
a longer time interval, requiring less frequent temporal 
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sampling. Another advantage of Ha over broad-band ob¬ 
servations, is that Ha images provide a greater contrast 
between novae and the bright galaxy background, result¬ 
ing in increased survey completeness, particularly near 
the bright nuclei of galaxies (e.g., Ciardullo et al. 1987). 

Based on the luminosity function of novae in both M31 
and the Galaxy, we estimate that approximately 50% of 
the Virgo novae should reach a peak absolute magnitude, 
Mmax,iia = —7.5®. Assuming a distance modulus to the 
Virgo cluster of /xq — 31.0 (see Table 1), we require ob¬ 
servations that can reach mua — 23.5. Numerical simu¬ 
lations revealed that with a mean seeing of 0.8" typical 
for the CFHT we could detect a m = 23.5 nova during 
grey time with a signal-to-noise ratio, S/N, ~ 12 at 60" 
from the center of its host galaxy and S/N ~ 6 at 30" 
from center in a one-hour exposure. To avoid saturation 
of the nucleus, and to enable the removal of cosmic rays 
hits, our 1 hr exposures were divided into 12 five-minute 
integrations. 

Our complete survey was spread over four epochs of ob¬ 
servations spanning two years. The first two epochs were 
separated by two months, which, given the persistence 
of the nova Ha emission, provides reasonable temporal 
coverage, with minimal overlap of detected novae. The 
third and fourth epochs were scheduled in the following 
year to ensure that we could discriminate between slowly 
evolving novae and other types of variable sources. 

Given an estimate for the absolute nova rate in M49 
(100 yr“^, Ferrarese et al. 2003), we conducted a simula¬ 
tion to estimate the number of novae we could expect to 
observe in each epoch. In the simulation, we consider de 
Vaucouleurs (1948) model for the brightness profile of the 
galaxy and then distribute artificial novae with a density 
that follows the background light. The peak brightnesses 
and fade rates for the artificial novae are based on Ha 
light curves of 14 novae observed in M31 (Shatter & Irby 
2001) and M81 (Neill & Shara 2004), adjusted to the 
distance of the Virgo cluster. Given the temporal sam¬ 
pling, exposure times (see Table 2) and an estimate of 
the nova rate we compute the number of novae expected 
to be observed above a given S/N. Under the null hy¬ 
pothesis that the nova rate is independent of GC specific 
frequency, and scales directly with K-band luminosity, 
we would expect to observe 22, 32 and 15 novae in M87, 
M49, and M84, respectively with a S/N > 10. On the 
other hand, in the extreme (and likely unrealistic) case 
that all nova progenitors are formed in GGs, the number 
of novae expected scales directly with the GG specific fre¬ 
quency and the predicted rates were 86, 32 and 7 novae 
for the three galaxies, respectively. 

3. DATA REDUCTION AND ANALYSIS 

In each observation, we positioned the target galaxy on 
the top of CCD chip 22 near the center of the MegaCam 
array. More than 99% of the total light from each galaxy 
falls on the same six central chips surrounding the nu¬ 
cleus (chips 12, 13, 14, 21, 22, 23), and these chips were 
the focus of our analysis. The images were not signifi¬ 
cantly dithered, and so chip gaps of ~ 13" persist in the 

® The Ho magnitude is defined on the AB system where 
ruHa = 0 for /;, = 2.53 X 10“® ergs cm“^ s“^ A“^. Through 
a 76A wide filter, this corresponds to a zero point flux of ^ 
1.9 X 10“^ ergs cm“^ s“l. 


final images for each epoch. The area obscured by chip 
gaps is ~ 3% of the total survey area. However, due to 
the positioning of the target galaxies, no chip gaps are 
encountered within a radius of 120" of the nuclei. Using 
the same photometry used to plot the cumulative back¬ 
ground light in Figure 2, we estimate the light enclosed 
within this radius to be greater than 60% of the total 
light for all targets, leaving considerably less than 1.2% 
of the light (and likely the novae) obscured by chip gaps. 
This loss is negligible compared to other uncertainties in¬ 
herent in nova rate calculations, and we have not made 
any formal correction for it in the analysis. 

The MegaCam data come pre-processed with bad pixel 
masking, bias subtraction, and flat fielding being per¬ 
formed prior to the data release. Once the pre-processed 
data were in hand, individual 5 minute exposures for each 
of the six central chips were spatially registered using the 
IRAF® routine WREGISTER. The individual CCD images 
were then median stacked using the tasks IMALIGN and 
IMCDMBINE to produce master images for each chip in 
each epoch, typically representing 60 min of total ex¬ 
posure time (see Table 2). Finally, the WCS for each 
master image was recalibrated to match between epochs 
using the IRAF task CCMAP in conjunction with precise 
coordinates for several reference stars on each image ob¬ 
tained from the US Naval Observatory (USNO) HI cat¬ 
alog (Monet 2003). 

3.1. Nova Detection 

To identify nova candidates the master images were 
spatially aligned, point-spread function (PSF) matched 
and differenced using the ISIS package (Alard & Lupton 
1998). In the resulting images objects of variable bright¬ 
ness are conspicuous, allowing us to detect them by eye. 
Variables qualified as nova candidates if they had an ap¬ 
proximately Gaussian PSF, and if they were visible in 
only one or two consecutive epochs, but did not appear 
on both sides of the 10 month gap in observations. 

Within the inner ~ 30" of each galaxy center the dif¬ 
ferencing became less effective due to the increasingly 
bright background light. In an attempt to detect nova 
candidates in the inner regions we applied the IRAF 
task MEDIAN, which produced a smoothed approximation 
of the background light profile. Stellar sources are not 
completely removed in the median smoothing process, 
but they are greatly reduced in amplitude. The result¬ 
ing median image was then subtracted from the original 
to produce a background-subtracted image of the cen¬ 
tral region where stellar sources are easier to identify. 
The background-subtracted images from separate epochs 
were then compared to identify nova candidates. No new 
nova candidates were identified using this method, how¬ 
ever we were able to confirm pre-existing candidates or in 
a couple of cases eliminate candidates due to non-stellar 
PSFs first identified in the background subtracted im¬ 
ages. All told, we have discovered a total of 27, 37, and 
19 novae in the fields of M87, M49, and M84. 

3.2. Photometry 

® IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatory, which is operated 
by Association of Universities for Research in Astronomy, Inc., un¬ 
der cooperative agreement with the National Science Foundation 
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To calibrate the nova Ha magnitudes, we began by 
performing aperture photometry using the IRAF task 
PHOT to determine instrumental magnitudes for the nova 
candidates and several nearby reference stars from the 
USNO B1 catalog, which were employed as secondary 
standards. For nova candidates close to the nucleus we 
needed to account for the steepness of the background 
light profile before we could perform accurate aperture 
photometry. To ensure that we were able to accurately 
measure the nova brightness, we pre-processed subimages 
centered on the nova using the task IMSURFIT, which ex¬ 
cludes a small circular region surrounding the nova and 
then calculates and subtracts a fit of the background 
light. 

Once instrumental magnitudes were in hand, stan¬ 
dard AB magnitudes for the novae (which correspond 
to monochromatic fluxes averaged over the Ha filter 
bandpass) were estimated through differential photom¬ 
etry with respect to the known i?-band magnitudes of 
the reference stars. Conversion to a true Ha magnitude 
that would reflect the total nova emission at Ha would 
require a knowledge of the width of the nova Ha emis¬ 
sion relative to the filter bandpass. Since the former 
value varies from nova to nova depending on the expan¬ 
sion velocity of the nova shell, an accurate measure of a 
nova’s Ha flux through filter photometry is not possible. 
For consistency with our earlier studies (e.g., see Ciar- 
dullo et al. 1990b; Shatter et al. 2000), we assume a 100% 
filling fraction in the present analysis. In this case, the 
adopted nova Ha magnitudes are simply equivalent to 
standard AB magnitudes, with ruHa = 0 corresponding 
to fx = 2.53 X 10“® ergs cm“^ s“^ A“^. 

Equatorial coordinates, the Ha magnitude, and the 
isophotal semi-major axis, a, are given in Table 3, 4, 
and 5 for the three galaxies, respectively. 

3.3. Spatial Distributions, Completeness, and Galaxy 
Membership 

The spatial positions of the novae detected in each of 
the three galaxies are shown in Figure 1. Given that the 
Virgo galaxies are members of a rich cluster, it is possible 
that some fraction of the novae, especially those seen at 
a large distance from their putative host galaxy, may in 
fact be associated with nearby dwarf galaxies, or possibly 
with intracluster novae (e.g., see Neill etal. 2005; Shara 
2006). To pursue this possibility further, and to explore 
how the nova populations are distributed in each of their 
host galaxies, in Figure 2 we have compared the radial 
cumulative distribution of the novae with the cumulative 
distribution of the galaxy’s background light using the g- 
band photometry of Cohen (1986) for M87 and M49, and 
the H-band photometry of Michard (1985) for M84. The 
different bandpass used for the M84 photometry should 
not affect the analysis significantly given that color gradi¬ 
ents in the Virgo ellipticals are not appreciable (e.g., see 
Cohen 1986). In each case the inner value of the isophotal 
semi-major axis, aim was set to be the minimum value at 
which novae could be detected. Artificial star tests dis¬ 
cussed below in section 4.2 show that our nova samples 
are complete to the limiting magnitude beyond Oin = V 
for M87, Oin = 1.4' for M49, and a™ = 0.5' for M84. 
Since our surveys were always incomplete within 0.5' of 
the center of each galaxy, we cannot address the possi¬ 
bility that the nova density is enhanced near the nuclei 


M87 (NGC4486) 



M49 (NGC4472) 



M84 (NGC4374) 



Figure 1. The spatial positions of the novae detected in each of 
our surveyed galaxies, plotted over equally spaced isophotes from 
/i = 20 to /i = 28 mag arcsec”^. Open circles represent novae that 
are more than 10^ from the nucleus, which have been excluded from 
the second analysis of M87. Filled red circles mark the locations 
of GC novae. 
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M87 (NGC4486) 



M49 (NGC4472) 



M84 (NGC4374) 



Figure 2. The cumulative distributions of the novae detected in 
each galaxy in our survey plotted against the cumulative distribu¬ 
tion of the background light. The inner value of a was set to assure 
the sample was was complete to the level of the faintest nova. The 
adopted inner radii of a = 1.0', a = 1.4', and a = 0.5', for M87, 
M49, and M84, resulted in the omission of 2, 7, and 4 novae from 
the full nova samples, respectively. KS statistics for the three fits 
indicate that the novae follow the background light in M87, M49 
and M84 with 33%, 66%, and 56% probability. 


of the galaxies, as was suggested by Madrid et al. (2007) 
in their HST study of novae in the inner regions of M87, 
and by Ferrarese et al. (2003) for M49. 

If the number density of nova progenitor binaries is 
simply proportional to the number density of stars, we 
would expect the cumulative nova distribution as a func¬ 
tion of isophotal semi-major axis to closely trace that 
of the galactic background light. Based on an examina¬ 
tion of the cumulative distribution plots, including their 
Kolmogorov-Smirnoff (KS) statistics, it appears that the 
density of novae follows the background light reasonably 
well. In the case of M87 there seems to possibly be two 
separate populations present. The M87 distribution is 
smooth from the central region to a ~ 10', then there 
is a gap in the distribution followed by what appears to 
be an increase in the nova density at a ^ 10'. A visual 
inspection of the M87 field reveals that within the region 
of this outer population several small neighboring galax¬ 
ies are being included in the survey, such as NGC 4478 
which by itself is ^ 10% the luminosity of M87. One 
of the novae in the outer population lies within ^ 2.7' 
of NGC 4478. There is also evidence for an extended 
halo surrounding M87 (e.g., see Kormendy et al. 2009). 
Thus is it possible that these outer novae are not directly 
associated with M87. 

Given the relatively poor fit of the nova distribution to 
the cumulative light in the case of M87 we have consid¬ 
ered a truncated sample of novae with a < 10' consisting 
of 21 novae for which there is a high probability of direct 
association with the galaxy. The revised cumulative dis¬ 
tribution is shown in Figure 3. After removing 6 novae 
with a > 10' and 2 novae with a < 1' and rescaling the 
remaining novae to the background light accordingly, the 
KS statistic between the background light and the nova 
distribution for M87 improves markedly (KS=81%). 

4. RESULTS AND DISCUSSION 

4.1. Relative Nova Rates 

As described earlier, the homogeneous nature of our 
survey, where we monitored each galaxy in a consistent 
fashion, allows us to estimate the relative nova rates sim¬ 
ply by comparing the relative numbers of novae detected. 
We have found that the relative number of novae dis¬ 
covered in each galaxy (27 in M87, 37 in M49, and 19 
in M84) is in good agreement with that predicted under 
the null hypothesis that the nova rate scales directly with 
galaxy luminosity. This is true whether we consider the 
entire population of novae discovered in the extended 
field surrounding each galaxy, or, in the case of M87, 
whether we consider the more restrictive sample (21 no¬ 
vae) based on the best match of the cumulative nova dis¬ 
tribution with the background light. Thus, based solely 
on the relative numbers of novae discovered in the three 
galaxies, we find no evidence that the GC specific fre¬ 
quency plays a significant role in the nova production in 
the host galaxies. 

4.2. Absolute Nova Rates 

Although our program was designed to compare rel¬ 
ative nova rates in the three Virgo ellipticals, it is also 
possible to use our data to estimate absolute nova rates 
in the galaxies as well. Our procedure follows that em¬ 
ployed in earlier studies where we used a Monte Garlo 
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Table 3 Table 4 

M87 Novae M49 Novae 


Nova 

N4486 

MJD 

a 

(J2000) 

S 

(J2000) 

"iHa 

(mag) 

a 

(') 

2011-01a 

55590.9 

12:30:26.4 

12:24:50 

22.8 

6.97 

2011-Olb 


12:30:36.9 

12:24:05 

22.2 

3.76 

2011-Olc 


12:30:43.5 

12:22:59 

22.0 

1.77 

2011-Old 


12:30:44.3 

12:24:53 

22.4 

1.91 

2011-Ole 


12:30:49.1 

12:23:08 

21.4 

0.35 

2011-Olf 


12:30:49.7 

12:22:04 

22.3 

1.44 

2011-Olg 


12:31:26.1 

12:21:06 

22.3 

11.12 

2011-03a 

55646.0 

12:30:23.8 

12:17:22 

21.2 

11.01 

2011-03b 


12:30:37.0 

12:21:47 

21.7 

4.49 

2011-03C 


12:30:46.0 

12:19:35 

22.4 

4.38 

2011-03d 


12:30:52.2 

12:22:13 

23.0 

1.42 

2011-03e 


12:30:52.3 

12:21:36 

22.3 

2.00 

2011-031 


12:30:56.3 

12:13:19 

22.7 

10.46 

2012-02a 

55983.0 

12:30:44.1 

12:24:01 

23.2 

1.50 

2012-02b 


12:30:44.4 

12:21:18 

22.1 

2.89 

2012-02C 


12:30:44.5 

12:23:03 

22.8 

1.46 

2012-02d 


12:30:47.3 

12:24:45 

22.0 

1.39 

2012-02e 


12:30:51.3 

12:23:47 

21.5 

0.59 

2012-021 


12:30:55.0 

12:20:55 

23.2 

2.90 

2012-02g 


12:30:56.2 

12:31:15 

22.9 

8.77 

2012-03a 

56006.9 

12:30:14.1 

12:35:50 

22.8 

15.33 

2012-03b 


12:30:28.9 

12:38:09 

22.4 

15.53 

2012-03C 


12:30:45.2 

12:23:54 

21.9 

1.18 

2012-03d=- 


12:30:51.2 

12:21:59 

22.5 

1.56 

2012-03e 


12:30:52.7 

12:25:06 

21.0 

2.03 

2012-031 


12:30:58.5 

12:22:57 

22.9 

2.60 

2012-03g 


12:31:26.0 

12:21:46 

22.7 

11.08 


Note. — The units for right ascension are hours, minutes, sec¬ 
onds. The units for declination are degrees, arcminutes, arcsec- 
onds. Apparent magnitudes are approximate only. The parameter 
a is the isophotal semi-major axis of the nova. 

This nova was detected coincident with GC M87: [P JB2009] 1667. 


M87 (NGC4486) 

Truncated at 10 arc min 



Figure 3. The cumulative distributions of the novae detected in 
our survey of M87 plotted against the cumulative distribution of 
the background light truncated to 10^ A KS test shows that the 
distributions match with 81% probability. 

approach to estimate the nova rates (e.g., Franck et al. 
2012). In this approach, for a given assumed intrinsic 
nova rate, R, we compute the number of novae that we 
can expect to observe, Vobs(.R), given assumed proper¬ 
ties of the nova population (peak brightnesses, and fade 
rates), our observing cadence, the distance to the galaxy 
in question and the limiting magnitude of our survey at 


Nova 

N4472 

MJD 

a 

(J2000.0) 

S 

(J2000.0) 

mua 

(mag) 

a 

(') 

2011-02a 

55593.9 

12:29:42.5 

07:57:58 

22.2 

2.63 

2011-02b 


12:29:43.4 

07:59:17 

21.9 

1.37 

2011-02C 


12:29:46.2 

08:01:25 

22.7 

1.44 

2011-02d 


12:29:47.7 

08:01:44 

22.9 

1.84 

2011-02e 


12:30:08.4 

07:58:17 

21.4 

6.50 

2011-03a 

55650.0 

12:29:30.7 

08:06:24 

21.5 

7.53 

2011-03b 


12:29:34.7 

07:57:30 

22.3 

4.87 

2011-03C 


12:29:37.7 

08:00:15 

21.2 

2.61 

2011-03d 


12:29:37.9 

07:59.04 

22.3 

2.91 

2011-03e 


12:29:42.1 

08:02:54 

21.4 

3.09 

2011-03f 


12:29:46.4 

08:04:41 

23.3 

4.87 

2011-03g 


12:29:47.2 

07:59:30 

21.7 

0.55 

2011-03h 


12:29:48.1 

07:59:32 

21.2 

0.62 

2011-03i 


12:29:50.3 

07:59:10 

22.9 

1.29 

2011-03j 


12:29:52.5 

07:58:58 

23.2 

1.89 

2011-03k 


12:29:52.7 

08:01:46 

23.2 

2.70 

2011-031 


12:29:53.2 

08:00:09 

22.6 

1.93 

2012-02a 

55983.1 

12:29:07.3 

08:13:10 

22.4 

16.57 

2012-02b 


12:29:36.6 

07:58:50 

21.9 

3.46 

2012-02C 


12:29:44.5 

07:59:21 

22.7 

1.01 

2012-02d 


12:29:47.4 

07:58:40 

22.2 

1.42 

2012-02e 


12:29:51.1 

08:00:04 

22.0 

1.33 

2012-02f 


12:29:56.9 

07:54:48 

22.5 

5.81 

2012-02g 


12:30:12.7 

07:55:59 

22.7 

8.26 

2012-03a 

56007.0 

12:29:18.4 

08:09:34 

23.5 

11.97 

2012-03b 


12:29:44.4 

08:02:36 

22.6 

2.66 

2012-03C 


12:29:44.7 

08:03:02 

22.5 

3.10 

2012-03d 


12:29:46.0 

08:00:45 

22.6 

0.74 

2012-03e 


12:29:47.9 

08:01:39 

22.0 

1.78 

2012-03f 


12:29:49.4 

08:02:08 

21.5 

2.42 

2012-03g 


12:29:50.9 

07:58:38 

22.6 

1.76 

2012-03h 


12:29:53.6 

08:00:51 

23.5 

2.28 

2012-031 


12:29:54.0 

08:00:23 

22.7 

2.20 

2012-03j 


12:29:57.0 

07:58:08 

21.8 

3.36 

2012-03k 


12:29:57.5 

07:56:46 

21.3 

4.30 

2012-031 


12:29:58.2 

07:58:23 

21.6 

3.58 

2012-03m 


12:30:06.1 

08:04:02 

21.9 

8.81 

Note. — 

- see notes for Table 3 




Table 5 

M84 Novae 

Nova 

N4374 

MJD 

a 

(J2000.0) 

5 

(J2000.0) 

(mag) 

a 

(') 

2011-02a 

55594.9 

12:24:40.8 

12:53:12 

23.4 

5.65 

2011-02b 


12:24:54.3 

12:52:38 

23.4 

2.43 

2011-02C 


12:25:07.1 

12:53:28 

21.8 

0.92 

2011-04a 

55652.8 

12:24:47.6 

12:52:07 

22.9 

4.19 

2011-046“* 


12:25:02.2 

12:53:43 

22.9 

0.63 

2011-04C 


12:25:03.5 

12:52:18 

23.0 

0.98 

2011-04d 


12:25:15.1 

12:52:28 

21.8 

2.88 

2012-02a 

55984.1 

12:25:00.3 

12:52:33 

23.7 

1.14 

2012-02b 


12:25:00.7 

12:53:49 

23.0 

0.95 

2012-02C 


12:25:01.0 

12:55:15 

23.2 

2.21 

2012-02d 


12:25:01.8 

12:53:48 

22.8 

0.75 

2012-02e 


12:25:02.6 

12:53:19 

22.3 

0.29 

2012-02f 


12:25:02.6 

12:53:32 

22.9 

0.42 

2012-02g 


12:25:19.6 

12:57:48 

22.7 

6.37 

2012-03a 

56007.9 

12:24:47.0 

12:58:38 

23.0 

6.87 

2012-03b 


12:25:04.2 

12:52:53 

21.9 

0.37 

2012-03C 


12:25:05.1 

12:53:22 

22.5 

0.42 

2012-03d 


12:25:06.2 

12:57:57 

22.0 

5.05 

2012-03e 


12:25:07.8 

12:53:09 

23.3 

1.03 


Note. — see notes for Table 3. 

^ This nova was detected coincident with a GC, NGC 
4374:[CLW2011] 033. 
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each epoch, i, as a function of position in the galaxy. For 
a range of plausible values of R, the simulation constructs 
a set of model Ha light curves by randomly selecting peak 
magnitudes and decay rates from a sample of actual Ha 
nova light curves from M31 (Shatter & Irby 2001) and 
M81 (Neill & Shara 2004). Given the dates of our ob¬ 
servations (see Table 2) and the distance modulus to the 
galaxy (see Table 1), an observed nova luminosity func¬ 
tion, ni(m,R), for each epoch is calculated. Following 
Franck et al. (2012), the number of novae we can expect 
to observe as predicted by our Monte Carlo analysis is 
then given by: 

-^MC(^) = EE C{m + Anii) ni{m, R), (1) 

i m 

where C{m) is the completeness as a function of apparent 
magnitude, which accounts for the variation in limiting 
magnitude across the galaxy. 

The completeness function is estimated through ar¬ 
tificial star simulations using the IRAF task ADDSTAR. 
Artificial novae are generated in several magnitude bins 
spanning the range of nova magnitudes observed. For 
each bin, the artihcial novae are distributed throughout 
the galaxy image, with a surface density that follows the 
background light. Then, employing the same search tech¬ 
nique used to discover the real novae, we determine the 
fraction of artificial novae that we can recover as a func¬ 
tion of magnitude. Completeness functions for the three 
galaxies are shown in Figure 4. Variation in the com¬ 
pleteness at different epochs, which may differ slightly in 
exposure time and seeing, is compensated for by estimat¬ 
ing an offset, Am^ = mumfi — between the epoch 

used in the artificial star tests and the other epochs. 

The most probable nova rate for a given galaxy is 
given when the number of novae actually observed, Vobsj 
matches the number of novae predicted by our Monte 
Carlo analysis, Nmc{R)- The Monte Carlo program re¬ 
peats this procedure 100,000 times and records the num¬ 
ber of times Nmc{R) = -^obs for each galaxy. The num¬ 
ber of matches as a function of R, normalized to the total 
number of trials, produces the probability distribution of 
nova rates for each galaxy shown in Figure 5. The peak 
of this distribution represents the most likely intrinsic 
nova rate, R, with Icr errors estimated by assuming Bi- 
Gaussian fits to the probability distributions (Buy & De 
Clerk 1972). We have also computed a nova rate for 
M87 based upon the truncated sample with a < 10'. A 
summary of the nova rates for each galaxy is presented 
in Table 6. We note that the rates determined for M87 
and M49, in particular, are somewhat higher that those 
found previously in studies by Shatter et al. (2000) and 
Ferrarese et al. (2003). Despite our higher nova rate for 
M87, our value of 154l^g yr“^ for the full sample falls 
short of the rate of ~ 200 — 300 yr“^ suggested by Shara 
& Zurek (2002) from early HST data, but is similar to 
that found by Mizusawa et al. (2013). 

In addition to computing nova rates, we have also de¬ 
termined the AT-band luminosity-specific nova rate, i^k 
for each galaxy from the data given in Table 1. In all 
cases we find values of vk of ~ 3 novae per year per 10^*^ 
solar luminosities in the K band. Our results are some¬ 
what higher, but still marginally consistent with the find¬ 
ings of Williams & Shatter (2004) and Giith et al. (2010) 
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Figure 4. The completeness functions, C{m), showing the frac¬ 
tion of artificial novae recovered in our artificial nova simulations 
for the different galaxies as a function of apparent magnitude. 














































































M87 (NGC4486) 



M49 (NGC4472) 



M84 (NGC4374) 



Figure 5. The probability distributions of the annual nova rates, 
R. The peak of each distribution represents the most likely nova 
rate, with la errors determined assuming that the distribution is 
Bi-Gaussian (see text). 

who found for a large sample of galaxies that i/x ^ 2 iLl 
independent of Hubble type. Our X-band luminosity- 
specific nova rates are presented in Table 6. 

5. GLOBULAR CLUSTER NOVAE 

Two of the novae discovered in our survey, one in 
M87 and one in M84, appear to be spatially coincident 


Table 6 

Nova Rates 


Galaxy 

.^^Obs 

R (pr-i) 

I'K 

(IQlO L0,x-J/r)-l 

M87 

27 

154123 

3.8 ± 1.0 


21 

124+2^ 

3.1 ± 0.8 

M49 

37 

189126 

3.4 ±0.6 

M84 

19 

95tl4 

3.0 ±0.6 


^ Only novae interior to an isophotal semimajor axis 
of 10^ have been considered in this analysis. 


with GCs associated with these galaxies. No GC no¬ 
vae were identified in M49. The GC nova candidates 
were discovered from the detection of residual light in 
differenced images at the positions of known GCs. In¬ 
strumental magnitudes were measured for the globulars 
in question, as well as for several nearby clusters, for 
all four epochs. While the magnitudes of the surround¬ 
ing clusters remained reasonably constant, the two GCs 
thought to host novae were seen to increase signihcantly 
in brightness in only one epoch. Figure 6 shows images of 
the M84 and M87 GC nova helds from epochs when the 
novae were detected, quiescent epochs, and their differ¬ 
ences. We attribute the residual light to erupting novae 
in these clusters. Nova magnitudes were estimated by 
subtracting the contribution of the host GC determined 
from quiescent epochs. The resulting nova magnitudes 
given in Tables 3 and 5, were consistent with those of 
the field novae, providing further evidence that the clus¬ 
ter brightenings were indeed the result of novae erupting 
from within the cluster. 

The PSFs of the GCs with erupting novae were consis¬ 
tent with those for other GCs, and were indicative of a 
single, unresolved light source. The centroid of the PSF 
for the GCs with erupting novae were coincident with 
their quiescent counterparts in all epochs to within 0.1", 
much less than the astrometric error of 1". We are suffi¬ 
ciently confident that these events are in fact GC novae 
to include them in our final count. 

The M84 GC, CLW2011:033, is a relatively blue clus¬ 
ter characterized hy g = 22.853 ± 0.016 and g — z = 
0.997±0.024 (Chies-Santos et al. 2011). The cluster ap¬ 
pears to be metal poor and relatively massive, being ~ 1 
mag brighter than the turnover magnitude. On the other 
hand the M87 cluster, PJB2009:1667, is characterized by 
V = 23.235T0.007 and V-I = 1.183T0.010 (Peng et al. 
2009). This cluster is redder than the M84 clster, more 
metal rich, with a flux that is ~ 0.5 mag brighter than 
the turnover magnitude. 

Despite predictions that novae may be common in 
GCs, relatively few have been reported. There have 
been just seven novae identified with GCs: four in M31 
[M31N 2007-06b in Bol 111 (Shatter & Quimby 2007), 
CXO J004345 in Bol 194 (Henze et al. 2009), M31N 2010- 
lOf in Bol 126 (Gao et al. 2012; Henze et al. 2013), and 
possibly an anonymous source in Bol 383 (Peacock et 
al. 2010a)], two in the Galaxy [T Sco (Nova Sco 1860) 
in M80 (Luther 1860; Pogson 1860), and possibly Nova 
Oph 1938 in M14 (Hogg 1964; Margon et al. 1991)], and 
an anonymous source in an M87 cluster (Shara et al. 
2004). Until now, the latter object was the only GC 
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Figure 6. Left column: Images of the M84 GC CLW2011-033 
from Epoch 2 (with nova, top), from Epoch 4 (without nova), and 
their difference. Right column: Images of the M87 GC PJB2009- 
1667 from Epoch 4 (with nova, top), from Epoch 1 (without nova), 
and their difference. In both cases residual light in the GCs is de¬ 
tected between epochs that we attribute to erupting novae. North 
is up and East to the left, with a scale of ~ 40” on a side. 


nova that has been reported in M87, with no cluster no¬ 
vae known in either M49 or M84. Clearly, the paucity of 
novae detected in GCs to date is in part a consequence 
of the difficulty in discovering novae against the bright 
background light of the cluster. 

The discovery of two GC novae from our survey pro¬ 
vides clear evidence that GCs do contribute some frac¬ 
tion of a galaxy’s nova system population, and that no¬ 
vae seem to be enhanced in GCs relative to the field. 
We can estimate the GC enhancement by estimating 
the mass in GCs in the three galaxies and comparing 
it with the integrated galaxy masses. The total number 
of GCs is estimated to be 13,200 in M87 (Harris et al. 
1998), 5900 in M49 (Rhode & Zepf 2004), and 1775 in 
M84 (Gomez & Rchtler 2004). Using a mean GG mass 
of 2.4 X 10® Mq appropriate for massive Virgo ellipti¬ 
cals (McLaughlin 1999), we estimate that the GG sys¬ 
tems have a combined mass of ^ 5 x 10® M©. For the 
galaxy masses, we adopt 7.6 x 10^^ M©, 8.1 x 10^^ Mq, 
and 4.0 x 10^^ Mq for M87, M49, and M84, respectively 
(Angeletti & Giannone 1997, Table 1), yielding a com¬ 
bined mass of ~ 2.0 x 10^® Mq. Thus, GGs represent 


^ 2.5 X 10“® of the total mass in the three galaxies. 
Given that we have detected two GC novae out of a to¬ 
tal of 83 novae discovered in our survey, we estimate that 
novae erupt ^ 10 times more frequently in GCs than they 
do in the field. 

Shara et al. (2004) estimated that novae occur in M87 
GCs with a frequency, / ~ 4 x 10“® novae per GC per 
year, which is ^ 2 orders of magnitude higher than what 
would be expected if novae were not enhanced relative to 
the field. In our observations of M87 we have discovered 
one GC nova out of 27 total novae. Assuming M87 has 
13,200 GCs and an overall nova rate of ^ 150 yr“^, we 
estimate a GC rate of / ~ 4 x 10““^ novae per cluster 
per year, or about an order of magnitude less than that 
estimated by Shara et al. (2004) for this galaxy. 

We can estimate the significance of the GC enhance¬ 
ment as follows. Let p = 2.5 x 10“® be the fraction of 
a galaxy’s mass in GCs. The probability of observing 
N or more GC novae out of a sample of M total novae 
assuming no GC enhancement is simply 


P- 


>N,M 


= E 

N=2 


Ml 


m{M -Ny. 




M-N 


( 2 ) 


The probability of observing 2 or more GC novae out 
of 83 novae is then P> 2,83 — 0.018, or ^ 1.8%. In other 
words, we conclude that novae are enhanced in GCs with 
98.2% confidence. Our estimate likely represents a lower 
limit to the true GC nova enhancement given the dif¬ 
ficulty in identifying novae against the GC background 
light. The selection effect is amplified when we consider 
that the most luminous GC have the densest cores, and 
are thus the most likely to form nova progenitor binaries. 

Novae discovered in GCs can provide insights into the 
dynamics of the clusters and isolate the effects of physical 
parameters such as age and metallicity in the (primarily) 
simple stellar populations of GCs. Low mass X-ray bi¬ 
naries, which are the only other binary systems that can 
be studied in Virgo cluster galaxies, provide an interest¬ 
ing point of reference. These compact neutron star and 
white dwarf systems are dynamically enhanced in GCs 
(Clark 1975; Hills 1976) by a factor of several hundred, 
and show several clear trends. The X-ray binary rate 
is correlated with the mass and metallicities of the host 
GCs. X-ray binaries are preferentially found the most 
massive GCs because the density of GCs increases with 
mass leading to a very strong correlation between the 
dynamical interaction rate, cluster luminosity and X-ray 
binary rate (Jordan et al. 2007; Peacock et al. 2010b). 
The fact that the two GCs that host novae in our survey 
are ^ 0.5 — 1 mag brighter than the turnover magnitude 
is consistent with this mass trend. We note, however, 
that because the peak luminosity of novae rival those of 
GCs it is very likely that some candidates in the most 
luminous GCs may be missed thereby underestimating 
the dynamical formation rate and the effect of GC mass 
(and the closely connected dynamical interaction rate) 
on nova formation. 

Extragalactic studies of X-ray binaries in GCs have 
firmly established that they are three times more abun¬ 
dant in the metal rich population as compared to the 
metal poor one in the typical bimodal GC populations 
in galaxies, a surprising effect that was hinted at in early 
Milky Way and M31 observations (Silk & Arons 1975; 
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Bellazzini et al. 1995; Kundu et al. 2002). The M84 GC 
nova found here is metal poor while the M87 candidate is 
metal rich. Interestingly, based on the published metal- 
licities or colors, of the 7 other likely or plausible nova 
GCs only two are likely to be metal rich. This may hint 
at the possibility that novae in GCs are not consistent 
with the metallicity effect seen in X-ray binaries. Clearly 
no strong conclusions can be drawn form this heteroge¬ 
neous set of observations with unknown biases. Future 
systematic observations such as this study are needed 
to test the intriguing possibility that more massive neu¬ 
tron star and black hole binaries respond differently with 
metallicity than WD binaries. 

6 . CONCLUSIONS 

We have reported the results of a Virgo cluster nova 
survey to test for a possible relationship between a 
galaxy’s nova rate and its GC specific frequency. In our 
CFHT survey spanning a total of 4 epochs over two years, 
we have discovered a total of 27, 37, and 19 novae in the 
fields of M87, M49, and M84, respectively. These data 
have led to nova rate estimates of of 154+^|, 189^22) 
95^14 per year, for the three galaxies. After considering 
the AT-band luminosities of the galaxies, we find similar 
luminosity-specific nova rates, i^k, for M87, M49, and 
M84 of 3.8 ± 1.0, 3.4 ± 0.6, and 3.0 ± 0.6 novae per year 
per 10^° Lq^k- These values are at the high end of the 
range of 1 — 3 novae per year per 10^° Lk ra found by 
Girth et al. (2010). 

The spatial distribution of the M87 novae shows some 
evidence that the nova population could perhaps be sep¬ 
arated into two distinct groups: that of an inner group 
which follows the background light out to a distance of 
about 10', similar to what is seen in M49 and M84, and 
an outer group of 6 novae clumped between ~ 14' and 
~ 16'. As mentioned previously, M87 appears to have 
an extended halo of diffuse light that is populated with 
lower luminosity elliptical galaxies such as NGC 4478. 
Thus it is possible that some or all of the novae seen in 
the outskirts of M87 may not in fact be associated with 
M87 itself. If we restrict our analysis of M87 to the 21 
novae with isophotal radii less than 10', we find that the 
cumulative distribution of novae closely follows the back¬ 
ground light distribution. For this truncated sample, we 
find a somewhat smaller nova rate of 124^47 per year for 
M87, corresponding to a luminosity-specific nova rate of 
vji = 3.1 ± 0.8, which is more in line with the values 
found for other galaxies. 

We have discovered two novae spatially coincident with 
known GCs: one in M87 and one in M84. During the 
course of our survey a total of 83 novae were discovered 
in M87, M49 and M84, suggesting that GC novae could 
comprise ~ 2.4% of the novae seen to erupt in these 
galaxies. This estimate is likely a lower limit given the 
strong observational selection against the discovery of 
novae in the brightest clusters (where they may be most 
likely to occur). Taken together, the GC systems com¬ 
prise only ~ 0.23% of the combined masses of the three 
galaxies. Thus, taken at face value, it appears that novae 
may be enhanced by at least an order of magnitude in GC 
environments. Additional monitoring of the GC systems 
of Virgo ellipticals, in particular M87, will be required 
before any definitive conclusions can be reached. 


In summary, the results of our survey are consistent 
with the hypothesis that the nova rates in the three Virgo 
galaxies studied here are to first order simply propor¬ 
tional to their mass in stars. Despite the detection of two 
novae associated with GCs in M87 and M84 suggesting 
that novae may be enhanced in GCs over the field, there 
is no compelling evidence that the overall nova rates in 
the Virgo ellipticals are particularly sensitive to the their 
GC specific frequencies, which vary significantly between 
the three galaxies. 
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